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Abstract 

The output spectrum of both gas and semiconductor lasers usually contains more 
than one frequency. Multimode operation in gas versus semiconductor lasers arises 
from different physics. In gas lasers, slow equilibration of the electron populations at 
different energies makes each frequency an independent single- mode laser. The slow 
electron diffusion in semiconductor lasers, combined with the spatially varying optical 
intensity patterns of the modes, makes each region of space an independent single-mode 
laser. We develop a rate equation model for the photon number in each mode which 
captures all these effects. Plotting the photon number versus pumping rate for the 
competing modes, in both subthreshold and above threshold operation, illustrates the 
changes in the laser output spectrum due to either slow equilibration or slow diffusion 
of electrons. 

1 Introduction 

The basic physics of mode competition in gas and semiconductor lasers has been known 
since the 1960's pQ- jl]. When the laser medium is in quasi-equilibrium and spatial variations 
in the mode patterns can be neglected, only a single frequency appears in the laser output. 
A rate equation model developed by Siegman [21 0] describes the subthreshold, transition, 
and lasing regions of operation for a single mode laser. Quantum theories of a single mode 
laser [Sj give similar (if not idential) results for the laser output power versus pumping rate 
as the simpler semiclassical models. 
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Figure 1: A four level lasing system representing both gas and semiconductor lasers. A 
current I pumps electrons from a lower thermodynamic reservoir (L) into an upper electron 
bath {U). Currents I a {Ib) then transfer electrons into the upper level of transition A 
(B). An electron electron scattering rate sab {sba) from the upper levels of transition A to 
transition B {B to A) is also present in this model. The upper lasing levels has occupation 
factor /a ifs), while the occupation factor of the lower lasing level is assumed to be = 
ih = 0). 
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In this paper we extend the rate equation models for a single model laser |21 El to include 
mode competition in both gas and semiconductor lasers. For the gas laser we develop a set 
of coupled rate equations which allow electron scattering between the lasing levels. The slow 
electron scattering rate in gas lasers allows the electron occupation factors in the different 
lasing levels to get out of equilibrium with each other, producing multiple frequencies in the 
laser output. In semiconductor lasers electron scattering between the different energy levels 
is rapid, keeping the occupation factors in the different lasing levels in equilibrium with 
each other. But spatial variation in the optical mode intensities inside the laser cavity favor 
different lasing modes in different regions of space within the semiconductor. Following 
Ref. |lj we extend the rate equations for a homogenous (semiconductor) line to allow for 
spatially varying mode patterns, generating multiple frequencies in the laser output. 

We assume both gas and semiconductor lasers are a four level system, with two intermedi- 
ate lasing levels, as shown in Fig. ^ For simplicity we consider only two lasing modes, mode 
A and mode B, with mode A the favored lasing mode. For the gas laser A and B represent 
localized atomic states, while for the semiconductor A and B are spatially extended states 
within the energy band of a quantum well. To simplify the mathematics we assume lower 
lasing level is always empty, having occupation factors /a = /b = 0. In the language of gas 
lasers this means we assume electrons in the lower lasing level empty very efficiently into an 
electron bath (r^ , — >• 0). In semiconductor language we would say there is extremely effi- 
cient hole capture into the quantum well. We therefore consider only pumping electrons into 
the upper lasing level. In physical systems would also have to guarentee charge neutrality 
while pumping the laser, and therefore also consider details of pumping and relaxation out 
of the lower lasing level. 

In this paper we limit consideration of laser mode competition to only 1-2 competing 
modes. In actual lasers many modes compete, and this case has been considered by Casper- 
son [U]. The rate equation models in this paper can be easily generalized to consider multiple 
lasing modes. 

2 Gas Lasers: Spectral Hole Burning 

We construct the rate equations for an inhomogeneous line following the example of a 
single mode laser from Refs. OEl- We consider a scattering rate sa^b = sab for electrons 
from mode A to mode B. If we firstly neglect optical transitions, the rate equation for the 
occupation probability /a for electrons in state A is 



Here Ia is the pumping current per state. Specializing to thermodynamic equilibrium (no 
pumping current) implies the rates sab and sba are related by a Boltzmann factor as 





Sab = Sba exp {Ea - EB)/kBT. 



(2) 
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The rate equation for the photon number ua in mode A is unchanged from that for a single 
mode laser 

^ = KA{nA + l)NA-^AnA. (3) 

Here Ka is an optical rate constant for the A transition, Na the number of A states (number 
of atoms of type A in a gas laser), and 7a the cavity escape rate for photons having frequency 
Ua- Putting Eqs. (jl}-® together (along with analogous equations for mode B) into a single 
matrix equation for the variables /a, fs, ^a, and gives 
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We now specialize to {Ea — Eb) ^ fc^T so that we can approximate sab — sba = s. 
As scattering rate s becomes large (s — > oo), the occupation factors /a and are forced to 
equal each other (/a = /b) in this approximation. Without the approximation sab — sba, 
and in the absence of any optical transitions, we would have the occupation factors forced 
towards a Fermi distribution f{E) having /a = f{EA) and = f{EB) as sab ^ oo. 
The approximation sab — Sba therefore makes only a minor correction to the occupation 
factors, and is therefore not essential for our analysis of mode competition. We therefore 
approximate sab — sba = s leading to 
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Inspection of the upper left quadrant of the matrix in Eq. (jSJ shows that the scattering 
rate s is negligible until s exceeds one of the spontaneous emission rates A a oi Ab- Thus, 
the transition from two independent lasers (s = 0) to a homogeneous line (s — * oo) occurs 
when the scattering rate s exceeds the spontaneous emission rates A a and Ab- This is true 
in open cavity lasers with luminescence through the sides of the laser cavity, so that the 
spontaneous emission rates Aa and Ab greatly exceed the cavity escape rates 7a and 7^- 
If the cavity is closed, so that no side luminescence occurs, the spontaneous emission rates 
are forced towards the cavity rates, i.e. Aa 7a and Ab 75. In the case of a closed 
cavity the transition from two independent lasers to a homogeneous line occurs when the 
scattering rate s is comparable to the cavity rates 7a and 7^. Note that on a homogeneous 
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line (s ^ oo) Eq. (jS)) simplifies to f = /a = fs and 
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where I = Ia + Ib and A = Aa + Ab- 

We solve Eq. (0) in steady state using an iterative technique. From the ith iteration for 
the variables f^, n\, and n^, we produce the {i + l)st iteration by 
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(7) 



We have tried several different types of initial guesses for fA-, Jb, ^a, and to start the 
iterative procedure, and the final results seem to be independent of the different initial 
guesses. The initial guess which seems to converge in the shortest time is to start in the 
subthreshold region and take for f\, n\, and the analytical results for two independent 
single mode lasers (s = 0). When incrementing to the next pumping rate, assume an initial 
guess for /a, /b, nA, and ub which are just the converged values at the previous pumping 
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Figure 2: (a) Photon numbers ua and ub and (b) normalized occupation factors fA/ftk 
Ib/ fth when the electron scattering rate between states A and B is s = 0. The iterative 
solution of Eq. (0) (solid lines) matches the analytical solutions for two independent single 
mode lasers (circles) having s = from Eqs. (jHI)-©- 

We use some results from the single mode laser [21 EI to frame our discussion of two 
coupled lasing modes. When the two laser modes are decoupled (s = 0), the threshold 
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Figure 3: (a) Photon numbers ua and ub and (b) normalized occupation factors fA/fth 
fs I fth when the electron scattering rate s between states A and B equals the spontaneous 
emission rate of mode A = \S)Aa). Scattering forces the occupation factors Ja and jn 
towards each other, reducing the threshold current for mode A and increasing the threshold 
current of mode B. 

currents for each mode are if^ = 'JaAa/NaKa and = 'JbAb/NbKb- The normalized 
pumping rate r is defined as r = Ia/II^- We define the ratio of the two threshold currents 
when the modes are decoupled as z = I^/I^. We assume the pumping current divides 
equally among the two states so that Ia = Ib = -^/2. The occupation factors above threshold 
are fixed at fA = fth = IaI^aKa and = /J^ = ■Jb/NbKb due to gain saturation. 

An important laser parameter is the number of luminescent modes p, where A a = pKa- 
For mode B this leaves the relation Ab = pzKb{'^a/1b){Nb/Na)- We further simplify by 
taking = Na = Nb and 7a = 7s = 7- In the absence of mode coupling (s = 0), the 
results for photon numbers versus normalized pumping rate r are then 



Figures EHU show the photon numbers and occupation factors versus normalized pumping 
rate r for different scattering rates s. We choose p = 10'' and 2; = 1.1 in Figs. I21E1 In 
Fig. |21there is no scattering between states A and B [s = 0), leaving two independent single 
mode lasers. Iterating Eq. ((7j) (solid lines) then just reproduce the analytical results from 
Eqs. ((HI)-© (circles) in Fig. |2fa). The normalized occupation factors fA/fth fs/ fth 
Fig. 121 increase approximately linearly with pumping below threshold and saturate above the 
lasing threshold. 

As the scattering rate increases to s = LOA^ in Fig. El and even further to s = S.OA^ 
in Fig. 15 we see the lasing threshold for mode A shifts to a lower pumping current. The 
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Figure 4: (a) Photon numbers ua and and (b) normalized occupation factors fA/ftt 
and fs/ fth when the electron scattering rate s between states A and B is five times the 
spontaneous emission rate of mode A {s = 5.0A^). The threshold current of mode A 
decreases slightly from Fig.|3l while the threshold current for mode B substantially increases. 

threshold current for mode B continues to increase as the scattering rate s increases. The 
increase in lasing threshold for mode B is much more pronounced than the decrease in 
threshold current for mode A. Inspection of the occupation factors for the two decoupled 
lasers in Fig. El^b) explains the threshold current shifts. Increasing the scattering rate s 
forces the two occupation factors towards each other. In Fig. El^b) we have fs > /a in the 
subthreshold region. Hence scattering between the modes will increase /a and decrease fs 
in the subthreshold region as seen in Fig. Efb) and Fig. El^b). Scattering then lowers the 
threshold current required for mode A to lase. Once mode A reaches the lasing threshold, 
additional scattering between states A and B makes it more difficult for mode B to raise its 
occupation factor to fs = required for mode B to lase. 

The threshold current for mode A can shift in either direction, up or down, with additional 
scattering between the modes. If fs > /a in the subthreshold region, the case we have 
chosen in Figs|2H additional scattering s lowers the threshold current for mode A. If the 
occupation factors obey Jb < /a in the subthreshold region, then the threshold current for 
mode A increases with additional scattering s. When the pumping current divides equally 
between the states A and B as we have assumed, the threshold current for mode A shifts 
down with additional scattering s if (f^/I^) > {fthl ^th)i equivalently, if the spontaneous 
rates obey Aa > Ab- Given our assumptions of 7^ = 7b = 7 and A''^ = Nb = N, we require 
Ka > zKb for additional scattering s to lower the threshold current of mode A. Since we 
choose z = 1.1 and Kb = 0.75Ka in Figs. |2Hthis condition is satisfied. Increasing z and/or 
Kb could reverse the inequality and raise the threshold current for mode A with increased 
scattering s. If the pumping current divides unequally as Ia = al and = (1 — a)I, 
the requirement for additional scattering s to lower the threshold current of mode A is 
(1 — a) A A > aAB- The threshold current required for mode B to lase will always increase 
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when we add additional scattering s (assuming mode A is the favored lasing mode with 

fth < fth)- 

In the hmit of s — >■ oo we move towards a homogeneous hne. Figure El shows the solution 
of Eq. (jS)) with with s = 100 A a- The photon number ub and occupation factor fs are 
now essentially fixed when mode A starts lasing due to gain saturation. Iteratively solving 
Eq. © produces essentially the same graph as shown in Fig. |S| The homogeneous line shown 
in Fig. El is the opposite limit of two independent laser lines shown in Fig. |21 Varying the 
scattering rate s interpolates smoothly between the solutions in Fig. |21 and Fig. |31 
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Figure 5: (a) Photon numbers ua and ub and (b) normalized occupation factors fA/ fth 
and fsl fth when the electron scattering rate s between states A and B approaches infinity 
(s = IOOAa). Since s ^ cxd we approach the limit of a homogeneous laser line and of a single 
mode laser. 
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3 Semiconductor Lasers: Spatial Hole Burning 



When there are spatial variations in the optical intensity in different modes, two lasing 
frequencies can coexist on a homogeneous line. Figure IHl shows the normalized optical mode 
intensities ImaP and I^bP for the lowest two longitudinal modes in a cavity. Near the center 
of the laser (region I), mode A is the favored lasing mode. However in region II, where the 
optical intensity > \ua\'^, mode B is the favored lasing mode. If the gain medium were 
confined to region I, only mode A would lase. Similarly, for the gain medium restricted to 
region II, only mode B would lase. For semiconductor lasers the gain media fills the entire 
laser cavity, so there is competition for the available optical gain between the lasing modes. 

Whether or not a single or multiple frequencies appear in the laser output spectrum 
depends on the size of the electron diffusion coefficient D. For single frequency laser operation 
to occur the electron must diffuse from region I to region II in Fig. IHl before the photon exits 
the cavity. If the photon escapes the laser cavity before the electron can diffuse from region 
I to region II, the regions are essentially independent as far as the laser light is concerned. 
Optically, the laser behaves as if two independent (single mode) lasers operate inside the 
cavity. The distance from region I to region II in Fig. IHl is approximately one quarter of 
the lasing wavelength (A/4). So for open cavity lasers we expect essentially single mode 
operation whenever D ^ (A/4)^A. If the cavity is closed (no side luminescence) so that 
A 'J A, the condition for single mode laser operation becomes D ^ (A/4)^7^. 




A=1 



B=2 



X 



I II 

Figure 6: Competition between two lasing modes A and B is possible on a homogeneous 
line when the optical mode intensities and |u_bP vary in space. Mode A is favored in 
region I, while in region II mode B is the favored lasing mode. 

To describe semiconductor lasers quantitatively we need to generalize Eq. (jH) to account 
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for the spatial variation in the mode patterns and electron density inside the laser. The 
occupation factors in each mode will be spatially varying such that /a('") and fs — > 

/^(r). If we introduce the position dependent density of states N{E,r), the total electron 
density is now 

p{r,t)=J2NiE,,r)fi{r,t). (10) 

i 

The electron density in states A and B are pA^r) = N{Ea, r)fA{r) and psir) = N{Eb, r)fB{r) 
The total number of active type A lasing levels is then 

Na= f N{EA,r)dV, (11) 

J a 

where denotes integration over that portion of the laser cavity containing the active lasing 
media. We further define the scattering rate per initial and final state density sab as 

= N{EA,t)N{Es,ry ^^^^ 

To account for spatial variations in the electromagnetic modes inside the laser cavity we 
introduce the mode functions UA{r) and ub{t) such that the electromagnetic energy density 
is given 

hujAnA{t)\uA{r)\^ = t{r)\EA{r,t)\\ (13) 

where EaIt, t) is the electric field of mode A and e(r) the dielectric constant. Since we must 
have 

J^e{r)\EA{r,t)\^dV = hujAriAit), (14) 

where the integration region L denotes the entire laser cavity, the mode functions UAir) are 
normalized as 

J^\uAir)\'dV = 1. (15) 

We can insert this factor of ' 1 ' from Eq. p5j) wherever necessary in order to generalize Eq. 
to account for spatially varying electromagnetic fields. 

Using Eqs. (fTUjl - ffT^j) . the generalization of Eq. (jT]) to account for spatial variations in the 
electron density and electromagnetic field intensity is 

^ = N{EA,r) {-pArsAB[N{EB,r) - ps] + PBSBA[NiEA,r) - pa]} 

+RAir, t) - [KAVL]nA\uA\^PA - AaPa + DaV^Pa- (16) 

Here RA{r,t) = N{EA,r)lA{r,t) is the total pumping rate per unit volume into the state 
A, Vl = Ji dV the volume of the laser cavity, and Da the diffusion constant of electrons in 
state A. The generalization of Eq. (jSI) to account for spatial variations inside the laser is 



duA 
"df 



[KAVL]{nA + l) I \uA\^pAdV --iAUA. (17) 
J a 
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Eqs. (fT^ -(fT7 j) can be used to construct a generalization of the coupled mode Eq. (jH) to 
account for spatial variations in the laser. 

Our interest is in semiconductors with homogeneous optical lines, so we do not pursue the 
full generalization of Eq. (jH). We assume the scattering rate sab —>■ oo in the semiconductor, 
so that we are back on a homogeneous optical line. We assume negligible separation of the 
energy levels as before so that {Ea — Eb) <^ ksT and s = sab = sba The occupation 
factors we therefore take to be in equilibrium with each other at each point in space so that 
/(r) = fA^r) = fB{r). With these assumptions we have 



+DV^p{r,t)+R{r,t). 



Here p = pa + Pb the total electron density, the total pumping rate is i? = Ra + Rb, 
and we have taken D = Da = -Db for the diffusion constant. The final coupled mode rate 
equations for a homogeneous semiconductor line that we solve are 



dt 

together with 
and 



^ = ~[kAVL]nA\uA?P - \kBVL]nB\uB?p - Ap + DV^p + R, (19) 



dUA 

"df 

duB 



{UA + 1) / [kAVL]\uA\^pdV - JAUA, (20) 

inB + l)JjkBVL]\uB\^pdV-^BnB. (21) 
The position dependent optical rate constants KA{r), KB^r), and A{r) are 

= ^^ MfJ.^\^^AlL V (22) 
N{EA,r) + N{EB,r) 

with 

N{EA,r} + N{EB,r} 

and 

A(r) = AAN(E„rn AsN(E, r) 

Eqs. (jl9|) -(|2H ) are similar to Eqs. (E.l.Qa) and (E.1.9b) for a single mode laser from 
Ref. P^. Eqs. (fT ^ - (pT|) should also be considered the generalization of Eq. (jH)) to account 
for spatial variations while lasing on a homogeneous line. We simplify further by taking 
the electron density of states to be constant in space so that the rate constants KA{r), 
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KB{r), and A{r) are independent of space. For simplicity and concreteness we consider 
competition between the longitudinal laser modes, though the same procedure would work 
for the inclusion of transverse cavity modes. We choose an Fabry-Perot type cavity having 
the normalized mode functions 

2 

|m^|2 = sm^{A7rx/L) (25) 

and 

= sm^{B7cx/L). (26) 
AlL 

Here Al is the cavity area, L the cavity length, A is the number of half wavelengths in mode 
A, and B the number of half wavelengths in the longitudinal cavity mode B. 



3.1 Slow Diffusion 

We solve Eqs. flT^ - (j7T|) in steady state using an iterative technique. We take the diffusion 
constant D = 0, letting us solve Eq. (fT^ for p and substitute back into Eqs. ((201)- (ED) leading 
to 
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In Eqs. (j2ZI)-(I2HI) we have used the ration of optical coupling constants z = Ka/Kb, the 
number of luminescent modes p = A/Ka, assumed equal cavity escape rates •ja = Ib, and 
defined the normalized pumping rate as r = R/Rf^ = RVa/lP- We now produce the (m+l)st 
iteration for the photon numbers from the mth iteration using 

/ , f {riAim) + l)r(x) sin^(A7rx/L)(t/x/Lj 

UAym + I) = p / ; r J— -— ; — ; r j- -— ; — -— , (29) 

JLa nA[m) sin [A-kx/ L) + [1 / zjUBym) sin [Btxx/L) + (p/2) 

and 

r {nB{m) + l)r{x){l/z) sin^Bnx/ L){dx/ Lg) 

' ^JL^nA{m)sm\A7Tx/L) + {l/z)nB{m)sm\B7rx/L) + {p/2y ^ ' 

Here Vg = A^Lg with A^ the cavity area and Lg the length of active media. Once we have 
iterated Eqs. (j^ - ljHOI) to convergence, we obtain the electron density from 

p + 2nAsm'^{A7rx/L) + {2/z)nBsm'^{B7rx/L)' ^ ' 

Here = ^Afthl^a = tZ-^aKi is the electron density when mode A reaches threshold in 
a single mode laser. For simplicity we also take the pumping rate r{x) to be a constant 
(independent of space). 
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Figure 7: Photon numbers ua and hb versus normalized pumping rate r when the optical 
mode intensities |maP and I^bP are (a) constant in space and (b) have spatial variation. 
Mode B cannot lase when the mode intensities are uniform in (a). Spatial variations in the 
mode intensities allow mode B to lase in (b). 

Figure [7| shows the photon numbers and ub for two modes competing on a homoge- 
neous line. In Fig. [7|^a) the optical mode intensities are constant, so that = = 
1/{AlL) (as we implicitly assumed for the gas laser of section |2)). Figure [7|^a) therefore mim- 
ics the case where spatial variations in the laser are negligible. Two other cases where we 
can neglect spatial variation of the mode intensities are in a ring laser or in a semiconductor 
laser with rapid electron diffusion. In Fig. d^a) the photon number in mode B is fixed 
whenever mode A begins lasing. The solution of Eqs. (j27|) - (j28j) therefore reproduces lasing 
on a homogeneous line whenever the optical mode intensities and {ubI"^ are constant. 
We let the optical mode intensities vary in space in Fig. [7|^b), where we have taken the 
lowest two longitudinal cavity modes {A = 1 and B = 2). Mode B can indeed begin lasing 
in Fig. [7|^b), but requires a higher pumping rate than for two independent lasers on the same 
optical line. We have chosen parameters 2; = 1.1 and p = 10^ in Fig. [7| The circles in Fig. [7| 
show the solutions from Eqs. (jHI)-© for two independent single mode lasers. 

Spatial vatiations in the optical mode intensities and |m_bP become less relevant 

when the optical rate constant Kb becomes small. The ratio of the rate constants z = 
Ka/Kb in Fig. [7|is z = 1.1. We increase z in Fig. |Hlto (a) z = 1.3 and (b) z = 1.5, raising 
the threshold current required for mode B to lase. For the parameter z = 1.5 in Fig. |Htb), 
mode B no longer lases for the range of pumping currents shown (0 < r < 3). Although 
spatial variations in the optical mode intensities are still present in Fig. |H1 they become less 
relevant when the optical coupling constant Kb for mode B is too weak. 

Figure IHl^a) shows the photon numbers ua and hb versus pumping and for two higher 
lying longitudinal modes having A = 6 and B = 7. The photon numbers ua and ub in 
Fig. El^a) are essentially unchanged from those for the two lowest cavity modes having A = 1 
and B = 2 in Fig. [7fb). Figure El and Fig. [7| use the same parameters, namely 2; = 1.1 
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Figure 8: Increasing the ratio of the optical rate constants from z = K^/Kb = 1-1 in 
Figure [7| to (a) z = 1.3 and (b) z = 1.5 increases the threshold current required for mode 
B to lase. Spatial variations in the optical intensities become less relevant when the optical 
rate constant for mode B becomes too small. 

and p = 10^. The photon numbers ha and ub versus pumping therefore have little (if 
any) dependence on the number of half wavelengths in the cavity. The weak dependence of 
Fig. Uta) on the number of half wavelengths (where A — B = 1) is because the fraction of the 
gain media where > [""Ap is essentially independent of the number of half wavelengths 
in the cavity, as can be checked numerically. 

Spatial holes are burned into the electron density p(x) in Fig. El^b), especially when mode 
B begins lasing. Figure IHfb) shows the electron density p{x) inside the active laser medium 
for different pumping rates r. Because the mode intensities have a node at the mirrors and 
there is no electron diffusion, the pumping rates can be read directly from the normalized 
density axis at the mirrors (points I = and I = 200) in Fig. E^b). The normalized pumping 
rates in Fig. [7|^b) are r = 0.6, 0.8, 1.05, 1.65, 2.25, 3.0. The electron density p{x) is essentially 
constant for pumping rates below threshold (r = 0.6,0.8) in Fig. [7|^b). There is a small 
variation in electron density for pumping rates below threshold, which is invisible on the scale 
in Fig. IHfb). Above threshold the variation in electron density becomes quite pronounced, 
especially when mode B begins lasing (r = 1.65,2.25,3.0). The growth of electron density 
p{x) as we move from the center of the gain media towards the mirrors is due to our neglect 
of diffusion. Since the optical mode intensites \ua\'^ and |u_bP have a node at the mirrors, a 
large spatial hole is also burned into the main body of the laser. Smaller spatial holes arising 
from the oscillating optical mode intensities produce oscillations in the electron density. 

3.2 Fast Diffusion 

When we include diffusion {D ^ 0), we can no longer solve Eq. p9|) directly for the 
density p{x). Instead we discretize the active lasing medium, taking lattice points xi = la. 
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Figure 9: (a) Photon numbers and (b) electron density p{x) for two higher lying longitudinal 
modes A = 6 and B = 7. The photon numbers versus pumping in (a) depend only weakly on 
the number of half wavelengths in the laser cavity. Spatial holes are burned into the electron 
density p{x) in (b), due both to the mirrors and the oscillating optical mode intensities. 



Here a is the lattice spacing and < / < Imax, with La = ImaxO' the length of the active 
medium. With this lattice Eq. (fT^ reads 
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Here t = D/a^ is the diffusion rate and 

-wi = KAVLnA\uA{xi)\^ + KBVLnB\uB{xi)\^ + A. (33) 

Given an initial guess for the photon numbers ha and n^, we can invert Eq. (p?^ for the 
density p{x) in steady state. Taking a hypothetical five point lattice we have 
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(34) 



We use zero derivitive boundary conditions to truncate the matrix in Eq. ()34|1 . After solving 
Eq. (jH^ numerically for pi, we can substitute this density back into Eqs. (| ^ - (PT|) to generate 
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updated photon numbers ua and in steady state. We move from the mth 
iteration for the photon numbers by 

n^(m + 1) = K(m) + 1) / \^S^\uA\^pdV, 

Ja 7^ 

and 

ns(m + 1) = (riBim) + 1) / \^^^\uB\^pdV. 

Ja 7b 
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Figure 10: (a) Photon numbers and (b) electron density p{x) when the electron diffusion 
constant is D = (9/400) A^A. Adding some electron diffusion has raised the threshold current 
for mode B and reduced spatial hole burning effects due to the mirrors. 

Figures ^1 and ^2 show the effects of adding electron diffusion to the photon numbers 
and electron density in Fig. El In Figs. ITUIITT] the number of half wavelengths in the cavity 
are A = 6 and B = 7 with a ratio of optical rate constants z = Ka/Kb = 1.1, the same 
parameters as in Fig. El Figure ^1 uses a diffusion constant D = lOOa'^A with L = 200a 
{t = lOOp with Imax = 200), while Fig. [TTl uses a larger diffusion constant D = 500a^A 
with L = 200a (t = 500j9 with Imax = 200). Since the cavity is six half wavelengths long 
(L = 6A/2), we have D = {9/A00)X'^A in Fig. UHland D = {A5/A00)X^A in Fig. im These 
values for the diffusion constant are in good agreement with our order of magnitude estimate 
for when diffusion should affect the laser output characteristics. 

Adding electron diffusion raises the threshold current required for mode B to lase, as 
can be seen in Figs. iTUT a) and fTTT a). The larger the diffusion constant, the greater is the 
threshold current required for mode B to begin lasing. Diffusion also reduces the spatial 
hole burning effects due to the mirrors, leaving only the smaller spatial holes due to the 
difference in the number of half wavelengths in the cavity between modes A and B. Some 
overall gradients are still visible in the electron density in Fig. ITUT b). while in Fig. fm b) the 
overall average electron density density is essential uniform (mirror effects are negligible). 



to the (m + l)st 
(35) 
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Figure 11: (a) Photon numbers and (b) electron density p{x) when the electron diffusion 
constant is D = (45/400) A^A. The diffusion constant is now large enough that mirror effects 
are negligible. The electron density oscillates essentially periodically and the photon number 
in mode B is essentially constant above threshold. 

Fig. HH^b) resembles the picture of electron density used to illustrate the effects of spatial 
hole burning in the laser in Ref. |2]. 
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4 Conclusions 



We have generalized the laser rate equations in Refs. Ej both electron scattering 
between the different lasing levels to describe spectral hole burning effects in gas lasers. In 
order to model spatial hole burning effects present in semiconductor lasers, and guided by 
Ref . jlj , we then further generalized the rate equation model to include the effects of spatially 
varying optical mode intensities in the laser. 

In order for multiple frequencies to lase simultaneously, either the energy spectrum or 
spatial variation of the optical gain must be broken up into many independent (single moded) 
lasers. Electron equilibration (scattering rate) is slow in gas lasers, and this allows the energy 
spectrum to be broken up into many independent frequency ranges. An order or magnitude 
estimate for single mode laser operation to occur in gas lasers is that the scattering rate 
between electrons in the different energy ranges must exceed the spontaneous emission rate 
(s ^ A). For semiconductor lasers the electron diffusion is slow, and the gain media can 
be viewed as many independent lasers at each point in space. Due to spatial variations in 
the optical mode intensities, different lasing modes will be favored at each point in space. 
Since the regions where different modes dominate lasing are spatially separated by about 
one quarter wavelength, we need the diffusion constant to exceed D ^ X'^A/16 for single 
moded operation in semiconductor lasers. Numerical simulations given in this paper agree 
with these two order of magnitude estimates for the transition from single to multiple moded 
laser operation. 

Finally, we can summarize some general (and well known) conclusions about single versus 
multiple moded laser operation. Firstly, all lasers are single moded for some range of pumping 
rates near threshold. The range of pumping rates for single moded operation is larger for 
more scattering between electronic states and for faster electronic diffusion. But a range 
of pumping rates for single moded operation nonetheless exists no matter how weak the 
equilibration or how slow the electronic diffusion (unless two degenerate states are lasing). 
Secondly, all lasers become multi-moded when pumped hard enough (unless the gain medium 
is first destroyed by too high of a pumping rate). Finally, bad economic analogies do not 
describe laser mode competition. Statements such as 'Laser mode competition is just like 
life. The rich get richer and the poor get poorer.' are clearly incorrect. Even as the photon 
number in mode A increases, the worst that can happen is that the photon number in mode 
B remains constant. Mode B can also begin lasing (become an economic success) either 
by electrons scattering from mode A (working in a supporting industry often created by a 
competitor) or by specializing its spatial mode pattern to take advantage of optical gain 
inaccessible to A (working in another area of the economy to exploit talents and resources 
unavailable to a competitor). 
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